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We show that matter effects change the — * Vt osciUation probabihty by as much as 70% for 
certain ranges of energies and pathlengths. Consequently, the i/^ —> survival probability also 
undergoes large changes. A proper understanding of i^^ survival rates must consider matter effects 
in P,iT as well as P^e • We comment on a) how these matter effects may be observed and the sign of 
A31 determined in atmospheric neutrino measurements and at neutrino factories and b) how they 
lead to heightened sensitivity for small &13. 

PACS numbers: 14.60.Pq,14.60.Lm,13.15.+g 



Two of the most outstanding problems in neutrino 
physics are the determination of the mixing angle ^13 |1| 
and the sign of the atmospheric neutrino mass-difference 
^31 Q' A knowledge of these parameters is crucial for 
understanding the form of the neutrino mass matrix. So 
far, most studies have concentrated on the — > v^. os- 
cillation probability P^c as the means of determining the 
above parameters j6| . This is because the passage of neu- 
trinos through earth matter dramatically changes P^c ■ 

In this letter we point out that the — > z/,- oscilla- 
tion probability P^^- can also undergo significant change 
(a reduction as high as ~ 70% or an increase of ~ 15%) 
compared to its vacuum values over an observably broad 
band in energies and baselines due to matter effects. This 
can also induce appreciable changes in the the muon neu- 
trino survival probability P^^ in matter. 

The muon survival rate is the primary observable in 
iron calorimeter detectors like MINOS 7] and the pro- 
posed MONOLITH 13 and INO @ and a major con- 
stituent of the signal in SuperKamiokande (SK) |^ , the 
planned BNL-HomeStake [l^l large water Cerenkov de- 
tector, and several detectors considered for future long 
baseline facilities. The r appearance rate as a signal for 
matter effects can also be searched for in special r de- 
tectors being thought of for neutrino factories We 
show that the energy ranges and baselines over which 
these effects occur are relevant for both atmospheric 
and beam source neutrinos for the above experiments. 
Since all matter effects sensitively depend on the sign of 
A31 and on 6'i3, observation of the effects discussed here 
would provide information on these important unknowns. 

Our discussion below uses the approximation of con- 
stant density and sets A21 = Agoi = 0. Consequently 
the mixing angle 9i2 and the CP phase 5 drop out of 
the oscillation probabilities. This simplifies the ana- 
lytical expressions and facilitates the qualitative discus- 
sion of matter effects. We have checked that this works 
well (upto a few percent) at the energies and length 
scales relevant here. However, all the plots presented 
in this letter are obtained by numerically solving the full 
three flavour neutrino propagation equation assuming the 
PREM 12 density profile for the earth. Further, these 



numerical calculations assume A21 — 8.2 x 10 ^ eV 



sin^6'i2 = 0.27 14J and 5 = |T1|. We consider mat- 
ter effects in neutrino probabilities only but discuss both 
the cases A31 = ±|A3i|. We find that dramatic matter 
effects occur only for A31 > 0. 

Review of P^o in matter: In vacuum, the v^. 



oscillation probability is 



P™^ = sin^ 023 sin^ 2013 sm^ (I.27A31L/E) , (1) 

where A31 = m| — m^ is expressed in eV'^ , L in Km and 
E in GeV. In the constant density approximation, matter 
effects can be taken into account by replacing A31 and 
013 in Eq. JQ) by their matter dependent values. 



A^"i = V(A3i cos 2013 - A)2 + (A31 sin 2013) 



^31 

sin 2013 = sin 2013 A31/A 



31 



(2) 



where A = 2-\/2GFnoE is the Wolfenstein term. The 
resonance condition is A = A31 cos 20i3, which gives 
Eres — A31 cos20i3/2A/2GF'T-e- Naivcly, one would ex- 
pect PJJJf* to be maximum at E = Eres since sin 205^ = 1. 
But this is not true in general because at this en- 
ergy A3" takes its minimum value of A3isin20i3 
and P™o''* remains small for pathlengths of L < 1000 
Km. P^o'^' is maximum when both sin20'i" — 1 and 
sin^ (1.27AfiL/E) = 1 = sin2[(2p + l)7r/2] are satisfied. 
This occurs when Ej. — r^mat 
[li: 



E™eak- "^^^^ gives the condition 



(2p-h l)7r5.18 X 10^ 
tan 2013 



Km gm/cc. (3) 



Here, p takes integer values. This condition is inde- 
pendent of A31 but depends sensitively on 0i3. Using 
the product pavL vs Lfor the earth (calculated using the 
PREM profile), where pav is the average density for a 
given baseline L, we identify the particular values of PavL 
which satisfy Eq. Q with p=0 for three different values 
of sin^ 2013. These occur at L~ 10200 Km, 7600 Km and 
11200 Km for sin^ 20i3 = 0.1, 0.2 and 0.05 respectively . 
Note that p=0 is the only relevant value of p in this case. 
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FIG. 1: P^e ,P/jT and P^^ plotted versus neutrino energy, E (m GeV) in matter and in vacuum for both signs o/ Am^i for two 
different baseline lengths, (a) for L= 9700 Km and (b) for L= 7000 Km. These plots use A31 = 0.002 eV^ and sin^ 26*13 = 0.1. 



given earth densities and baselines p^. 



Matter effects in P 



In vacuum we have 



P™; = cos4(?i3sin2 2^23sin2(1.27A3iL/E), 
= cos^ 6*13 sin^ 26*23 sin^ (I.27A31L/E) 



(4) 



Including the matter effects [T^ changes this to 

P^f = cos^ 6l'°sin2 26*23 sin^ [1.27(A3i + A + A;J1)L/2E] 
+ sin^ 6li'^sin2 26*23 sin^ [1.27(A3i + A - A;j°i)L/2E] 

(5) 
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Compared to P™^', the matter dependent mass eigen- 
states here have a more complicated dependance on the 
Vfj, and Vr flavour content. Labeling the vacuum mass 
eigenstates as vi, 1^2 and i>3, v\ can be chosen to be al- 
most entirely and 1^2 to have no v^^ component [19| . 
Inclusion of the matter term A leaves V2 untouched but 
gives a non-zero matter dependent mass to v\. As the 
energy increases, the v^. component of i/™ decreases and 
the , Vt components increase such that at resonance 
energy they are 50%. Similarly, increasing energy in- 
creases the i/q component of (and reduces z^^, Vt com- 
ponents) so that at resonance it becomes 50%. Thus in 



the resonance region, all three matter dependent mass- 
eigenstates z^f , and contain significant and Vt 
components. 

We seek ranges of energy and pathlengths for which 
there are large matter effects in P^t, i-e for which 



AP„ 



"pmat 



T3vac 



is large. We show that this occurs 



for two different sets of conditions, leading in one case to 
a decrease from a vacuum maximum and in another to 
an increase over a broad range of energies, 
(i) Large decrease in P™.^* in the resonance re- 
gion : At energies appreciably below resonance, the 
cos2 6'5^ term in Eq. © « P™'^ (since 6*5^:^ = 6*13, 
A << A31, Ag^ K, A31) and the sin^ 6*5^ term is nearly 
zero. As we increase the energy and approach reso- 
nance, cos^ 6*5^ begins to decrease sharply, while sin^ 
increases rapidly. However, if resonance is in the vicinity 
of a vacuum peak, then the decrease in the cos^ term 
has a much stronger impact on P™^* than the increase 
in the sin^ 6*5^ term, since the latter starts out at zero 
while the former is initially close to its peak value (~ 1). 
As a result, P™^* falls sharply. This fall is enhanced by 
the third term in Eq. (Q, which is essentially 0.5 x PJfj** 
(which is large due to resonance ), leading to a large over- 
all drop in P™r* from its vacuum value. Note that the 
requirement that we be at a vacuum peak to begin with 
forces AP^T- to be large and negative , with the contribu- 



tions from the first and the third term reinforcing each 
other. 

The criterion for maximal matter effect, Eios — Ep^^jj., 
leads to the following condition: 

pV^r"" - (2p+ l)7r5.18 X 10^ (cos 2^13) Km gm/cc. 

(6) 

Unlike Eq. ((Hjl, which has a tan 2^13 in its denominator, 
Eq. lO has a much weaker dependence on 9i^. This 
enables one to go to a higher value of p without exceeding 
the baselines relevant for observing earth matter effects. 
Incorporating the Eics = Epcak condition we get AP^t- as 



cos 



sin 26*13 (2p 



1)- 



- 1 



(7) 



where we approximated cos 2^13 ~ 1. We note that, in 
general, AP^,- will be larger for higher values of both 
p and 013. From Eq. lO, for p=l and sin^ 2^13 = 

9700 Km ( 9300 
For 



peak occurs at 
and AP^^ w -0.7 (from Eq. 0) 



0.1(0.2,0.05), E,es =E 
Km, 9900 Km 
p=0, Eq. ® gives UJ^^"" - 4400 Km for sin^ 26*13 = 0.1. 
However, AP^t- is roughly one-tenth of the p=l case. In 
general, for a given baseline, the choice of an optimal p 
is also dictated by the constraint that the vacuum peak 
near resonance have a breadth which makes the effect 
observationally viable. 

In Fig. n^a) we show all three matter and vacuum 
probabilities for 9700 Km. In these plots A31 is taken as 
0.002 eV2, which gives E^os = E™^^^ - 5 GeV. The mid- 
dle panel of Fig. Q^a) shows that near this energy P™^* 
(~ 0.33) is appreciably lower compared to P™'^ {'-^ 1). 
Thus the drop due to matter effect is 0.67, which agrees 
well with that obtained earlier using the approximate ex- 
pression Eq. {T)). 

In Fig. 13 we show the 6I13 sensitivity of F'^f at 9700 
Km. In particular, at Ercs — E^^j^ the strong depen- 
dance on 6*13 is governed by Eq. {T)) above. Unlike PJ^Jf*; 
where the event-rate decreases as 6*^3 for small values of 
013, the T appearance rate at Ej-es ~ E™^]^ increases with 
decreasing 6*13. As sin^ 26*13 goes from 0.2 to 0, P™^* 
varies from 0.05 to ~ 1. For very small values of 
sin^20i3(< 0.05) it will be impossible to see a maximal 
resonance enhancement in P^^ because the distance for 
which this occurs exceeds the diameter of the earth. How- 
ever, the observation of resonant suppression in P^^. is 
possible, even for very small values of ^13, if the criterion, 

Wr(013 = 0) - iV,(0i3) > 3(VA^,(0i3 = 0) + ViV,(0i3)), 

is satisfied for the tau event rate. 

In general the resonance has a width, and this fact af- 
fects observability. To include the width of the resonance, 
we write A = A31 (cos 26*13 + q sin 2^13). We find that the 
large matter effects discussed above still do occur as long 
as A is within the width of the resonance or — 1 < g < 1. 



(ii) Increase in P^^ away from resonance 



It 



is also possible for P™.^ to differ appreciably from P™'^ 
away from resonance. This is evident in Fig. ^a) (cen- 
tral panel) in the energy range 7.5 GeV - 15 GeV. This 
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FIG. 2: P^T- plotted against the neutrino energy, E(in GeV) 
for different values of On ■ We have used A31 = 0.002 eV^. 



effect is an enhancement rather than a drop, i.e., AP^^- 
is now positive. AP^o is small in most of the latter part 
of the energy region under consideration and does not 
contribute in an important way overall. The dominant 
contribution to this enhancement arises from the sin^ ^"3 
term in P^,- (Eq. (O) which is large for i? >> Ercs- Since 
(A31 -I- A — A™j) « 2 A31 for these energies, we obtain a 
enhancement(~ 15%) which follows the vacuum curve. 
The difference between the vacuum and matter curves 
largely reflects the difference between the cos^ 6113 multi- 
plicative term in the vacuum expression Eq. and the 
sin^ ^13 multiplicative term in Eq. (0 • While this effect is 
smaller compared to the effect in (i) above, it occurs over 
a broad energy band and may manifest itself in energy 
integrated event rates. 

Finally, we comment on the observability of the matter 
effects in P^t- The energies in question are above, but 
close to the r production threshold. This suppresses the 
T appearance rates, and will necessitate a high luminosity 
beam experiment. Such direct observation must perhaps 
await the advent of superbeams and/or neutrino facto- 
ries. However, the effects in P™.^* manifest themselves 
indirectly in P™^' , as we discuss below, and these can be 
observed in an atmospheric neutrino experiment. 
Matter effects in P^^ : The deviation of P™^' from 

P™'^ clearly results from the combined effects in P™^* 



and P-<f' I.e. AP^^ 



AP, 



-AP - 

In case (i) above, for instance, AP^,- is large and negative 
while AP^o is positive and hence they do not contribute 
in consonance. However, the resulting change in P^,^ is 
still large, given the magnitude of the change (« 70%) in 
P^T- . This is visible in the bottom panel of Fig. Q^a), in 
the energy range 4-6 GeV. 

One also expects a significant drop in P™^' when either 
of AP^^r or AP^c is large and the other one is small. The 
first of these cases (AP^,- large, AP^e small) is shown in 
Fig. n^a) in the energy range 6 — 15 GeV, with the 
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enhancement in P™^' reflected in the decrease in P™^'. 
The second case (small AP^^, large AP^o) occurs when 
a minimum in the vacuum value of P^,- resides in the 
proximity of a resonance, and even the rapid changes 
in sin^ and cos^ 9fi^ in this region fail to modify this 
small value significantly. This condition can be expressed 
as I.27A31L/E = pTT. Note that this corresponds to a 
vacuum peak of P;^^- Substituting E as E^cs gives the 
distance for maximum matter effect in P^^ as 

pL'^i^'' ~ p TT X lO'^ (cos 2013) Km gm/cc (8) 

For p=l this turns out to be ^ 7000 Km. This effect 
[2lj | is shown in the bottom panel of Fig. Qlb). The large 
(40% at its peak) drop in P^^ seen in this figure derives 
its strength from the resonant enhancement in P^c ■ A 
sensitivity to 613 around the peak similar to the one dis- 
cussed above for P™^* also exists here, leading to a larger 
muon survival rate as 613 becomes smaller. 
The width of both these effects is significant, ranging 
from 4—10 GeV in the first case (Fig. ^b)) and 6—15 
GeV in the second. We have checked that they persist 
over a range of baselines (6000 - 9700 Km), making them 
observationally feasible. 

Observational Possibilities and Conclusions: We 

have shown that large matter effects in neutrino oscilla- 
tions are not necessarily confined to t'^ —^ Vc or Vj- 
conversions, but can be searched for in — > oscilla- 
tion and — > survival probabilities. We have dis- 
cussed their origin by studying the inter-relations of all 
the three matter probabihties, P™''*j ^^t'' and P™^*, and 
identified baseline and energy ranges where they act co- 



herently to give observationally large effects. The effects 
discussed are strongly sensitive to the sign of Am^^ , as is 
apparent in the figures above. Also, there is sensitivity 
to small 6*13 at the energy and baseline ranges identified 
for Pj;V''* and P;^^*. 

Specialized r detectors operating in long basehne sce- 
narios ^3 should be able to observe effects like the ones 
discussed in the central panel of Fig. Qfa) and in Fig. |2] 
Similarly, detectors capable of measuring muon survival 
rates, e.g. magnetized iron calorimeters can detect the 
effects visible in the bottom panels of Fig. Q^a) and (b) 
|20l |. To illustrate the observability of the effect, we calcu- 
late that for a magnetized iron calorimeter detector 0| 
and an exposure of 1000 kT-yr, in the energy range 5-10 
GeV and Lrange of 6000 - 9700 Km, with A31 = -1-0.002 
eV^ and sin^ 26*13 = 0.1, the total number of atmospheric 
/i~ events in the case of vacuum oscillations is 261. How- 
ever, it reduces to 204 with matter effects. The rates for 
/i+ in matter are identical to the vacuum value of 105 
events. These numbers reflect a Aa signal for the effect 
discussed above for P™^* [l^ . The Fermilab to Kamioka 
proposal has a baseline of 9300 Km and is within 
the range of baselines where these effects are large and 
observable. 

Finally, we remark that although the effects discussed 
here appear only for v^^ for A31 positive (and only for 
for A31 negative) it may still be possible to search for 
them in the accumulated SK data. 
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